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The spectrum of the integrable spin-S Heisenberg chains is completely characterized in terms of 
spin-T? spinons. In the continuum limit they form a quasi-particle basis to the higher level SU(2) 
Wess-Zumino-Witten (WZW) models. Enumerating the spinon states in finite systems we obtain 
effective single particle distribution functions for these objects which generalize Haldane's generalized 
exclusion principle to quasi-particles with non-Abelian exchange statistics. 
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Haldane's notion of fractional exclusion statistics jl|] has been used successfully to characterize elementary excitations 
in low dimensional quantum systems. For the excitations in the lowest Landau level arising in the fractional quantum 
Hall effect (FQHE) this concept has been shown to coincide with the usual definition of fractional statistics, namely 
through their properties under exchange of anyons g]. A complete description in terms of quasi particles obeying 
fractional statistics is possible for integrable quantum systems where detailed knowledge of the spectrum allow for an 
classification of all states in terms of the elementary excitations. Their statistical properties are described in terms of 
the statistical interaction parameters g a p which determines the many particle Hilbert space dimension 
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Here 

d« = dW- £ 9ap(N -6 aP ) (2) 


is the dimension of the Hilbert space of a single particle of species a in the presence of N a — 1 other particles of 
that species at fixed coordinates, are constants which can be interpreted as the dimensions of the single particle 
Hilbert spaces. This definition includes bosons (g a p = 0) and fermions (g a p — S a p). An new feature in this picture are 
the possible off diagonal statistics parameters g a ^p encoding the mutual statistical interaction of different species of 
anyons. This is realized in spin-i Heisenberg chains |lj and has recently been discussed in the context of the statistics 
of the quasiparticles and -holes in the FQHE at filling v = 1/m with m an odd integer |}|4|] . 

The approach outlined above allows to introduce single particle distribution functions which already contain all 
information on the statistical properties of the elementary excitations of the system thereby giving an easy access 
to thermodynamic properties of the system such as the low temperature specific heat. Unfortunately, it is easily 
seen that Haldane statistics (||) fails to describe correctly the interesting case of quasiparticles obeying non-Abelian 
statistics as realized for example in the FQHE at half integer fillings : the distribution functions obtained within 
this picture always lead to a low temperature specific heat characteristic of a Gaussian conformal field theory (CFT) 
with central charge c = 1 while the quasi particles for the v = i FQHE are described by a SU(2) 2 WZW model 
with c = |. 

Recently, an extension of Haldane statistics to quasi particles obeying non-Abelian exchange statistics has been 
proposed based on recursion relations for truncated conformal spectra constructing the basis of the CFT 

explicitely in terms of quasi-particle operators simple distribution functions could be calculated that coincide with 
the ones found from (||) for Abelian statistics and correctly reproduce the central charge for the non-Abelian case. 

In this letter we apply this method to characterize the spinon excitations of the Bethe Ansatz soluble spin-S* 
SU(2) Heisenberg chains |p^ , |TT| . The particular form of the Hamiltonian is not of interest here, it contains exchange 
interactions of nearest neighbour spins only that are given as polynomials of degree 2S in Sj ■ Sj+i. The eigenstates 
for the spin S chain of length L are characterized by complex numbers \j , j = 1, . . . , M satisfying the system of Bethe 
Ansatz equations 
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It is well known that regular Bethe states (corresponding to solutions with Xj < oo) are highest weight states with 
spin Stot = LS — M, limiting M/L by S [^2|. Moreover, in the thermodynamic limit L — > oo all solutions of (||) are 
built from so called strings A^") + \ (n + 1 — 2j) with real centers A^ n ^ and j = 1, . 



,71. 

Based on this classification the completeness of Bethe states has been proven for the periodic chain after 
rewriting (||) in terms of the real A^") any Bethe state is characterized uniquely by a collection of integer or half 
integer numbers ■ Now the number of possible states built from v n n-strings is just the number of admissible 

choices of quantum numbers ^ for j ^ k. The latter are found to be restricted to the intervals [— Qmax, Qmix] 
where 



oo 

E 
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2Q^ ax + 1 = L mm(n, 2S) + v n - Y] G mn v m (4) 



with G mn = 2min(m,n). In particular, the antiferromagnetic vacuum of a chain of even length is the unique singlet 
state in which all L/2 modes for the 25-strings are occupied. 

In a recent study Jl4|] of a related ideal system the "strings" appearing as solutions of a simplified version of the 
Bethe Ansatz equations (^) have been interpreted as the elementary objects spanning the Fock space of the system. 
In this picture one obtains the infinite dimensional matrix G mn for the statistical interaction of the system (see also 
p5|). One should note, however, that the role of the various strings in building the excitations over the physical 
vacuum, namely the antiferromagnetic ground state of the system is quite different: in the thermodynamic limit, 
the energy and momentum of a configuration is completely determined by the distribution of holes in the sea of 25* 
strings, while the ones with n ^ 2S are merely responsible for the correct counting of states fliof . In particular, 
different configurations of the n > 2»S-strings produce the 2 m -fold degeneracy of the multiplet due to the spin ^ 
degree of freedom carried by the m holes (spinons) while the n < 2S strings distinguish states with the same number 
of holes and total spin. Their presence is reflected in the S'-dependence of the critical behaviour of the system jll] : 
the low-energy effective field theory of the system has been identified as a SU (2)2S WZW model. Separating this 
field theory into a bosonic sector (with central charge c — 1) and a Z^s parafermionic sector []l6| the n < 2S strings 
have been assigned to the latter. In a different approach they have been identified as spanning the state space of a 
restricted solid on solid (RSOS) model on the basis of the number of m-particle states Jl7| . 

First let us briefly recall the situation in the S = | case where the ground state is a completely filled band of 
v\ = L/2 (real) 1-strings. The elementary excitations above the antiferromagnetic vacuum can be identified with 
holes in their distribution ]l2] | . For chains of even (odd) length L states with even (odd) number m of holes form 
"Yangian" multiplets containing states with total spin S — O(^), . . . ,m/2 which differ in the number of n > 1-strings 
fixing the total spin of the resulting state. Counting these states (including their SU(2) multiplicity) the dimension 
of the m-spinon subspace is 

d m =( L + 1 )- (5) 

Summing over the admissible m this gives the correct total number of states 2 L . Here the statistics of the excitations is 
modeled by hard core lattice bosons reflecting the level-1 SU (2) Kac- Moody symmetry of the critical theory. Taking 
into account the spin degree of freedom one finds that the spectrum of the m-spinon states decomposes exactly into 
the states present in the tensor product of m spins ^: in the m = 2 spinon sector there are ( L ^ +1 ) triplets (built from 
1-strings only) and 2 ) singlets (containing a single 2-string) . Thus we are led to distinguish spinons with different 
spin projection and to rewrite (^|) as a sum over configurations with nf (n^ = m — n\) spinons with spin f (J.) as 
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dm ^ +7i| ,m jj^ 
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where d a (n^ njj = (g) (L — — nj,) + 1. This gives g aa i — ^ establishing the semionic nature of the spinons just 
as in the ideal spinon gas realized in the Haldane-Shastry model p],p^|. The difference between this model and the 
nearest neighbour Heisenberg model considered here is that a spinon-spinon interaction lifts the degeneracy of states 
with different total spin in the super multiplets when their density becomes finite. Following Rcf. Jl9| the distribution 
function of spinons is 

2 
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which gives the low temperature specific heat C/L — itT/3vf of a c = 1 model (vf is the Fermi velocity of the quasi 
particles) . 

For S > i one finds from ([I]) that there is a one to one mapping between the possible configurations of n > 2S 
strings between the spin-5 model and the spin-i Heisenberg chain (see Table [|): for any m-spinon Bethe state with 



total spin < m/2 in the S = | chain with string configuration {vn^ 2 ^} there exist Bethe states in the spin-S' chain with 



string configuration such that f^s+n— 1 = v ™ ^ or n > 1 with the same number m of holes in the distribution 

of 2S , -strings. In addition, there are certain n < 2S strings present in these states. Note that this is seen only when 
considering holes, the number v 2 g of 2«S , -strings differs for different configurations of the shorter strings. Proceeding 
as in the S — h case the number of 2m spinon states in a chain of even length L is found to be 
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with no = m, n 2 s = in the product. Eq. (Q) and a similar expression for odd L correctly reproduce the total number 
of states for the spin-S' chain. However, an interpretation in terms of Haldane statistics, i.e. the ansatz ^ for spin-| 
quasi particles with an additional internal degree of freedom turns out to be impossible. 

To proceed we apply the method introduced in ||,^(| . We identify the Bethe states of the spin chain with the spinon 
basis of the chiral SU (2)k =2 s WZW model proposed in (a realization of the chiral CFT is the open boundary 
version of the Takhtajan-Babujian models whose eigenstates are enumerated by the same relation (0)). In this basis 
the m-spinon states are of the form 
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where each operator <^*a. . creates a spinon with spin projection er^ =y, J, and energy e(rii + Ai) = (2ttvf / L)(jii + Ai) 
( 1 0) is the Fock vacuum). The numbers ji in (H) take values 0, |, . . . , | subject to the fusion rules \]\ — \ = 5 of 
the CFT and we have defined = A(j^) — A(j„-i) with A(j) = j(j + l)/(fc + 2). A basis of the spinon Fock space 
is provided by the states (||) with o\ = . . . = cr^ =|, (Jn^+i = ■ ■ ■ = ca^+a^ —i with the modes ni = n^ m i n + hi 
and non decreasing sequences of non negative integers {hi}^2 1 and {^ij^^+V The minimal allowed mode sequence 
is constructed from ni. m i n = and 
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"i, min 



+ 1 ifjj+i =ji-i <3% 
otherwise 



To compute the distribution function for these modes we introduce truncated partition functions q A ^ X^\q, x, h) 
of the chiral CFT for which the occupied spinon modes should satisfy < £ and = j |8|,^(| . Here /1 + ah are 
chemical potentials for the spinons and we write q = e~P^ VF / L \ x = e^. These partition functions satisfy recursion 
relations 
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with initial condition X, 
given in terms of D(y) = diag(l — y, . 



(J) 




CjX 3 , Cj 



sinh((j + l)/3h/2)/smh{(]h/2). The matrix TZe = R(xq e , (3h) D(x 2 q 2e - 1 ) is 
, 1 — y, 1) and the symmetric matrix 



R(y,0h) 



/I 



V 



ciy 



C2V 



(l+c 2 y 2 ) (c iy + c 3 y 3 ) 

(1 + c 2 y 2 + c 4 y 4 ) 



C2SV 2S 

c2s-iy 2S ~ 1 + c 2S +iy 2S+1 
c2s-2y 2S ~ 2 + c 2S y 2S + c 2S+2 y 2S+2 



v^2S 2 
Ei=0 C 2iV 



J 



(10) 



Below it will be useful to factorize the matrix TZe in © further as lZelZ i _i with TZe = R{xq e , (3h)D{x 2 q 21 1 ). The 
non-zero matrix elements of this new matrix are (R(y, /3/z.))^ fe . = Ci-jy l ~i for < i — j < k. 
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Remarkably, the number of m-spinon states contributing to the partition function X\ ' coincides with (Q) for the 
spin-5 chain of even length L = 21, Since the finite size spectrum of low lying excitations of the spin chain is 
determined by the operator content of the SU{2)k WZW model we can identify Xi ' with the partition function 

Zl of the spin chain of length L — 2£ and similarly q^ k ^xi k ^ = Z 2 t-i- In this sense the partition function of the 
integrable Heisenberg chain can be generated by a recursion relation similar to (^) with the matrix TZe- In addition, 
this observation provides further evidence for the equivalence of the spinons of the integrable chain with the ones 
forming the quasi particle basis (||) to higher level SU(2) WZW models of Ref. [pl| . 

Viewing the k + 1 spinon modes added in the £ th iteration of (|J) as a single degenerate level in the single particle 
spectrum we can now approximate the partition function of the spin chain as Zn — Yli=i Af with the maximal 
eigenvalues A-j = A^ + '(y = xq l ,(3h) of the recursion matrix R(y, f3h)D(y 2 ). This leads to a spinon distribution 
function 

n(e) = 2yd y In A<+> (y, (3h) , y = e' 13 ^ (11) 

for the spin-5 chain. For k = 1 this function has been studied in Ref. [^) where complete agreement with the 
corresponding results from Haldane exclusion statistics (e.g. (||) for h = 0) could be established. In the general 
non-Abelian case k = 2S > 1 it is not possible to obtain a closed expression for n(e), for large negative e we find 
n(e) ~ 2k — (2 + 4<5/ Ci 2)e 2/3 ' e_ ' i * ) . The maximal occupation 4S of a level containing 2S + 1 spinon modes coincides with 
the result obtained from ([?]) using the fusion rules of the CFT. Furthermore we have checked numerically for spins up 
to S — 3 that the distribution function (|ll|) gives the correct behaviour JllJ] of the low temperature specific heat [^2| . 

For k — 2 the eigenvalues determining the distribution function for h = and h — * oo (leading to complete 
polarization of the spinons) are the largest real solutions of 

(A-l) 2 =y 2 (A + l) for/i = 0, 

(A-1) 2 (A + 1) =y 2 e^/ 2 A 2 for h - oo. (12) 

In Fig. [l] the distribution functions obtained from (|l^) are presented in comparison with the one for particles obeying 
g = i exclusion statistics. 

To summarize, considering integrable lattice realizations of the SU{2)k WZW models we have shown that the 
truncation scheme for chiral spectra of a CFT introduced in B can be interpreted as a recursion relation for partition 
functions of the corresponding lattice model. The truncated partition functions of the field theory show complete 
equivalence between the well established spinon picture for the integrable spin chain and the spinon formulation of 
the CFT. For the spin chains this method allows to compute distribution functions for spinons obeying non-Abelian 
exchange statistics. Frequently, the low-lying excitations in one-dimensional S — \ systems have been interpreted by 
invoking a spinon-picture (see e.g. Ref. p3|). Similarly, we expect that the characterization of the spinons in higher 
S models given here will prove useful for a better understanding of massive perturbations of the latter such as the 
S = 1 Haldane system. 

The authors thank K. Schoutens and J. Fuchs for discussions. This work has been supported in parts by the 
Deutsche Forschungsgemeinschaft under Grant No. Fr 737/2-2. 
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FIG. 1. Distribution function obtained from (|ll|) for the spinons of the spin-1 chain: shown are n(e) for vanishing magnetic 
field (bold curve) and for the fully polarized spinon system (dashed curve) in comparison with the result for g = -x exclusion 
statistics (thin curve). 



TABLE I. String configurations for a Bethe state with given number m of holes in the distribution of 25* strings and total 
spin Stot: [r^p 2 ] denotes a solution with v n = 1, v p = 2. For S — 2 the configuration of strings with n > 2S (< 2S) are listed 
separately, for a Bethe state the n > 2S configuration has to be supplemented by one of the possible n < 2S configurations 
listed in the last column. 



m Stot 


s=k 


S = 2 


2 1 



[2 1 ] 


[5 1 ] [3 1 ] 


4 2 
1 



[2 1 ] 


mm 

[5 1 ] [3%[2 1 ] 
[6%[5 2 ] [S 2 ],^] 


6 3 
2 
1 



[2 1 ] 

[4 1 ],t3 1 2 1 ],t2 3 ] 


[5 1 ] [s'us^'ui 1 ] 

[n^ 1 ]^ 3 ] [3 3 ],t3 1 2 1 ],tl 1 ] 
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